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ABSTRACT 

The  stability  properties  of  six  different  adaptive  schemes  with 
respect  to  model  order  error  are  analyzed.  Bounds  on  parameter  identi¬ 
fication  and  state  errors  are  established.  All  adaptive  schemes  considered 
are  robust  in  the  sense  that  the  error  is  of  order  of  the  "speed  ratio  u" 
between  the  modeled  slow  phenomena  vs.  the  neglected  fact.  The  dependence 
of  the  error  on  the  input  signal  is  shown  to  be  crucial.  The  bounds 
obtained  indicate  possibilities  for  reducing  the  error  by  a  proper  choice 
of  the  input  signal. 
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1.  INTRODUCTION 

Global  stability  properties  of  model  reference  adaptive  systems  [1]- 
[10]  are  guaranteed  under  the  "matching  assumption"  that  the  model  order  is  not 
lower  than  the  order  of  the  unknown  plant.  Since  this  restrictive  assumption  is 
likely  to  be  violated  in  applications,  it  is  important  to  determine  the  robust¬ 
ness  of  adaptive  schemes  with  respect  to  such  modeling  errors.  Recently  several 
attempts  have  been  made  to  formulate  and  analyze  reduced  order  adaptive  identi¬ 
fiers  [11].  The  results  of  such  studies  depend  on  the  characterization  of  the 
model-plant  mismatch . 

In  this  paper  we  examine  stability  properties  and  performance  of 
various  types  of  identifiers  and  adaptive  observers  [1]-[10]  when  the  model- 
plant  mismatch  is  due  to  a  fast  ("parasitic")  part  of  the  plant,  and  the  order 
of  the  model  is  equal  to  the  order  of  the  slow  ("dominant")  part  of  the  un¬ 
known  plane .  We  express  our  results  in  terms  of  a  "speed  ratio"  u  of  the  slow 
versus  the  fast  phenomena.  Scalar  y  is  small  and  positive  and  u  i  0  means  that 
the  fast  part  of  the  plant  reaches  its  steady-state  instantaneously,  that  is 
the  plant  order  reduces  to  that  of  its  slow  part.  The  fictitious  "reduced- 
order"  plant  it  thus  obtained  when  in  the  actual  plant  y  >  0  is  replaced  by 
y  *  0. 

This  singular  perturbation  approach  is  a  convenient  parameterization 
of  the  model-plant  mismatch.  In  our  formulation  adaptive  observers  are  designed 
for  the  reduced  order  plants,  but  they  are  applied  to  the  actual  plants.  In 
Section  2  we  derive  a  singularly  perturbed  state  space  realization  of  the  plant 
and  give  a  statement  of  the  problem.  In  Section  3  we  analyze  the  stability 
properties  of  an  identification  scheme  [1],  [2]  employing  a  lower  order  model 


» 


and  we  obtain  bounds  for  parameter  identification  and  state  errors.  In  Section 
4  we  analyze  the  stability  properties  of  reduced  order  minimal  adaptive  observers 
[5]- [7]  applied  to  actual  plants  and  obtain  bounds  for  parameter  and  observation 
errors.  The  stability  properties  of  a  nonminimal  adaptive  observer  [4],  [8] 
designed  for  a  lower  order  plant  and  applied  to  the  actual  higher  order  plant  are 
analyzed  (Section  5)  and  bounds  are  obtained  on  the  parameter  identification  and 
output  errors.  In  Section  6  similar  results  are  obtained  for  the  parametrized 
adaptive  observer  [9].  A  qualitative  analysis  based  on  these  bounds  is  given  in 
Section  7  illustrated  with  computer  simulation  results.  Particularly  important  is 
the  sensitivity  of  the  parameter  identification  er-or  with  respect  to  the  excita¬ 
tion  input  signal. 

2 .  PROBLEM  STATEMENT 

Systems  possessing  .>* .  and  fast  parts  can  be  represented  in  the 
explicit  singular  perturbation  form 


x  »  Anx 

+  A12xf  +  Bxu 

(2.1) 

uxf  =  A21x 

+  A 22xf  +  B2u 

(2.2) 

y  “  c 'x 

(2.3) 

where  x,  x^  are  n  and  m  vectors  respectively,  u  is  an  r  control  vector  and  u  is 
a  small  positive  parameter  associated  with  the  presence  of  "parasitic"  elements, 
such  as  time  constants,  masses,  etc.  [14].  The  matrices  A^.A^  .A^^.A^^  ,B^  and  B2 
have  appropriate  dimensions. 

Without  altering  the  input-output  characteristics  of  the  system  we 
will  use  the  transformation  [13]  7,  *  x^  +  Lx  +  A^B^u  and  analyze  the  equivalent 


representation 


kz 


3 


D 


x  ■  Ax  +  Bu  +  HT] 

(2.4) 

■V 

ur|  -  +  uA^Bg  u 

(2.5) 

« 

t- 

y  *  c'x 

(2.6) 

* 

where 

A  •  A^  -  A^L,  A^  *  ^2  ^  uLA^2»  ^ 

-  B2  +  uLBj 

— 

L  satisfies  the  algebraic  equation 

A22L  -  A21  +  uLA^2L  ”  mLah  *  0 

(2.7) 

* 

Approximate  expressions  for  L,  A,  A^  and  B  are 

.  -1 

• 

L  *  A22A21  +  0(u) 

(2.8) 

■ 

A  *  A^^  -  A^2A22  .+  0(u)  ,  Aj  •  A22  +  0(li)  , 

(2.9) 

Bj  -  b2  +  oo*),  b  -  bl -a12a‘^b2  +0(t*) 

(2.10) 

r 

► 

Representation  (2,4)  -(2.6)  containing  is  found  to  be 

% 

convenient  for  getting 

tighter  error  bounds  and  clarifying  the  dependence  of 

the  error  on  the  characteris 

tics  of  the  input. 


i 


Dominant  Part 


i 


.  1.  Representation  of  the  plant 


The  part  of  the  system  described  by  (2.4)  and  (2.6)  will  be  referred 
to  as  the  dominant  part  of  the  plant  whereas  (2.5)  will  be  called  the  parasitic 
part  of  the  plant.  Suppose  that  an  adaptive  scheme  is  designed  for  the  n-order 
dominant  part  assuming  that  there  are  no  parasitics  ie  HTj ■  0.  This  scheme  is 
then  applied  to  the  actual  plant  with  parasitics.  The  purpose  of  this  paper  is 
to  examine  the  robustness  of  the  scheme  with  respect  to  the  parasitic  part  of 
the  plant  and  to  obtain  bounds  on  the  parameter  and  output  or  3tate  errors. 
Throughout  the  paper  the  following  assumptions  are  made: 

(i)  A  is  stable  and  A^  is  asymptotically  stable 

(ii)  The  order  of  the  dominant  part  of  the  plant  is  known 

(iii)  The  tripple  (A,B,C)  is  completely  controllable  and  completely 
observable 

(v)  The  only  available  signals  are  u(t)  and  y(t) 

(vi)  u(t)  and  u(t)  are  piecewise  continuous  bounded  functions 
of  time. 

It  will  be  shown  in  the  following  sections  that  the  stability  of 
several  adaptive  algorithms  in  the  presence  of  parasitics  is  equivalent  to  the 
stability  of  a  linear  time-varying  equation  with  a  parasitic  input 

Z(t)  -  A  (t)Z(t)  +  Hll(t)  (2.11) 

n 

where  Z(t)  is  a  composite  error  vector.  It  should  be  pointed  out  that  (2.11)  is 

not  input  n  to  state  Z(t)  linear  because  A  (t)  depends  on  x  which  in  turn  de- 

n 

pends  on  n.  This  dependence  will  be  explored  for  each  particular  scheme.  Our 
approach  is  to  first  derive  conditions  under  which  the  homogeneous  part  of 
(2.11)  is  uniformly  asymptotically  stable  (u.a.s.)  for  each  n  of  interest. 


5 


After  these  conditions  are  found.  Lemma  1  is  used  to  obtain  bounds  on  Z(t). 
Lemma  1:  If  the  homogeneous  part  of  (2.11)  is  u.a.s.  then  Z(t)  is  bounded. 
A  bound  on  the  norm  of  Z(t)  as  t -►«  is  of  order  of  u  and  is  given  by 


lim  »Z(t)ll  <uy  —  —  Hh#  B A~1B_fl 
m2  °2  f  f 


(2.12) 


Proof :  Since  (2.11)  is  u.a.s.  there  exist  positive  numbers  m^  and  m2  such  that 
its  transition  matrix  <j>(t,t)  satisfies 

-m.(t-t) 

H  <p(t,t)B  <m^e  for  all  t>x  and  all  x>0 


Therefore  from  (2.11)  we  can  write 


t  -m2(t-x) 


HZ(t)U<m^e  II Z ( 0 ) U  +  /gin^e  U Hll  fi  n ( t) H dx 


(2.13) 


Since  A^  is  asymptotically  stable  and  u(t)  is  bounded  by  assumption  we  set  from 


(2.5) 


-V/„  t  -a,^ 

H  n  ( t )  n  <  a-  e  B  n  (0)  u  +  f  a.e  II  Ap  BJ  ydx 

1  q  1  i  r 


(2.1'; 


where  y  *  sup  Ou(x)H  and  a.,  a.  are  positive  constants.  From  (2.13)  and  (2.14) 
t>r  1  1 


we  have 


I Z ( t) 8  <  uy  ^  ^  IHHA^BJ  +  m.e  Z(0) »  -  a  m 

m2  a2  f  f  1  1  1  (ij-  a2) 


a.  m.  .  a1  II APB  J  H  Hll  a.,  ra  H  HW  -a0t/u 

UY  ~~  ~~  I  Hi  I  A,  B  J  +  UY  -7-  m  —  - }+  ■  —■  e  2  {3  n  (0>| 

a2  m2  f  f  “2  1  ,  °2,  ,  a2, 

(m2  -  — )  («2  -  — ) 

II  A^BJ 

"  WY  - 

? 


(2.15) 


and  (2.12)  follows  as  t-*-». 
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Bound  (2.12)  is  convenient  because  it  w'  1  be  shown  that  factor 
ml  °1  -1 

uy - I  A,  B 9  H#  remains  the  same  for  all  adaptive  schemes  considered  in  this 

m2  <*2  f  f 

paper.  The  dependence  of  (2.12)  on  confirms  that  the  schemes  are  robust  because 

a 

Z^O  as  The  factor  ~  IIa^  B^IHIhII  is  determined  by  the  parasitics,  m^ 

depends  on  the  initial  error  Z(0)  and  depends  on  the  rate  of  convergence  with¬ 
out  modeling  error.  The  presence  of  parameter  Y  characterizing  the  input  will  be 
shown  to  have  a  crucial  effect  on  the  parameter  and  state  errors.  In  Sections 
3,  4,  5  and  6  we  establish  u.a.s.  of  the  homogeneous  part  of  (2.11)  and  derive 
specific  forms  of  the  bound  (2.12)  for  six  different  adaptive  schemes. 


3.  IDENTIFICATION  UU2] 

It  is  desired  to  identify  the  pair  (A,B)  in  (2.4)  by  using  an  nth 
order  model  and  assuming  that  the  state  vector  x  is  available  for  measurement. 
The  presence  of  the  parasitic  input  HTj  is  disregarded  in  the  design  of  the 
identification  algorithm. 


The  nth  order  model  for  the  identification  of  the  pair  (A,B)  is  given 


by  [1],  [2] 


x  ■  K  (x  -x)  +  A  (t)x  +  B  (t)u 
mm  mm 


(3.1) 


where  K  is  a  stable  matrix  and  the  adaptive  laws  for  adjusting  A  (t)  and  B  (t) 

m  tn 


-iiex' 


(3.2) 


f  *  -r2eu' 


(3.3) 


where  $  »  A  (t)  -  A,  Y  £  B  (t)  -  B  and  e  >  x  -  x  are  the  parameter  and  state 
mm  m 


errors  and  >0,  I"2  -  ~'2  >0. 


J 


I 


1 1 


In  the  absence  of  parasitics  (q«0)  it  is  shown  in  [2]  that  if  u(t) 
is  sufficiently  rich  for  an  nth  order  plant  (ie  the  components  of  u(t)  are 
linearly  independent  and  each  component  contains  at  least  distinct  frequen¬ 
cies)  then  e,  <p,  V  +  0  as  t-*-00.  The  stability  of  the  identification  algorithm 
in  the  presence  of  parasitics  is  equivalent  to  the  stability  of  the  following 
system. 


e  -  Ke  +  <|>x  +  4fu  -  Hq 


(3.5) 


<!>  *  -T^ex' 


(3.6) 


4*  =  — T^eu ' 


(3.7) 


where  (3.5)  is  obtained  by  subtracting  (2.4)  from  (3.1).  To  express  (3.5)-(3.7) 

as  a  linear  time  varying  equation  in  the  form  of  (2.11)  we  define  Z(t)  *  [ef  ,<p’  ,f’  ]  ' 

where  $  =  ,<J>  ]’,  ?  =  ]’  and  <fr.,  V .  are  the  ith  rows  of 

l  l  n  1/  n  li 

p  and  ¥  respectively.  Then  we  denote 


An(t) 


K  I  x’  0  ;  u\.  0 
i  0  **x*  !  0  -u* 

-r.  • 

lx  ■ 


-»1 


(3.8) 


where 


Y(1)X 
Y1  1 


Y2  X2 


Y(1)x 
n  n 


Yi  ui 


Y2  u2 


(2) 

Yr  ur 


and  y^^  is  the  jth  row  of  i*l,2.  To  apply  lemma  1  we  now  investigate 


whether 
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Z(t)  -  An(t)Z(t) 


(3.9) 


is  u.a.s.  or  not.  The  stability  of  (3.9)  can  be  easily  established  by  choosing 
the  same  Lyapunov  function  as  in  [1],  [2]  for  the  case  without  parasitics  (n  •  0) . 
However  for  u.a.s.  the  components  of  the  vector  have  to  be  linearly  indepen¬ 

dent  functions  of  time.  A  sufficient  condition  for  the  case  n  •  0  is  that  the 
components  of  u(t)  be  linearly  independent  and  each  component  contains  at  least 
distinct  frequencies.  In  our  case  x  depends  on  n  which  in  turn  depends  on 
u.  Thus  in  some  cases  n  might  destroy  the  "richness"  property  of  u.  This  can 
be  avoided  by  choosing  u  such  that  the  components  of  the  vector  [x'.u'.n']'  are 
linearly  independent  functions  of  time,  for  which  a  sufficient  condition  is  that 


the  components  of  [u'.u']'  be  linearly  independent  and  each  component  of  u  be  the 

n+m+1 

sum  of  sinusoids  with  at  least  — j —  distinct  frequencies.'  This  implies  that 
the  components  of  [x',u']  are  linearly  independent  functions  of  time  and  (3.9)  is 
u.a.s.  Thus,  lemma  1  immediately  furnishes  bound  (3.4),  since  HhO  *lH#.  We 
summarize  this  result  in  theorem  1. 

Theorem  1:  If  u  is  sufficiently  rich  for  the  (n-Hn)-th  order  plant  and  the  compo¬ 
nents  of  [u'.u']'  are  linearly  independent  then  the  identification  algorithm 
(3.1)-(3.3)  is  stable  in  the  sense  that  the  composite  error  Z(t)  is  bounded.  The 
bound  is  of  order  of  y  and  is  given  by  (2.12). 

Remark  1;  To  guarantee  the  u.a.s.  of  (3.9)  it  is  sufficient  to  make  u  sufficiently 
rich  for  the  highest  suspected  order  of  the  actual  plant.  Although  this  is  a 
feasible  approach  in  most  applications,  there  is  a  considerable  "overkill"  in 
requiring  this  richness.  In  fact  it  can  be  shown  that  the  system  (3.9)  will 
remain  u.a.s.  for  almost  all  u  which  are  sufficiently  rich  for  the  nth  order 
dominant  part  of  the  plant  only.  For  example  each  component  of  u  can  contain  any 
distinct  frequencies  except  for  a  particular  combination  for  which  the  condi¬ 


tion  of  linear  independence  of  x  and  u  can  be  lost. 


.  MINIMAL  FORM  ADAPTIVE  OBSERVERS  [7 


The  plant  (2.4)-(2.6)  Is  assumed  to  be  single  Input  single  output  and 
an  nth  order  adaptive  observer  Is  designed  to  estimate  the  state  vector  x  of 
the  dominant  part  of  the  plant  and  to  Identify  the  triple  (A,B,C)  or  Its  equiva¬ 
lent.  The  presence  of  the  parasitic  input  H**!  Is  disregarded  during  the  design. 
The  stability  of  the  adaptive  observer  operating  on  the  real  plant  in  the  pres¬ 
ence  of  the  parasitic  input  is  then  analyzed.  Two  different  types  of  minimal 
adaptive  observers  are  considered  separately.  Case  1  and  Case  2. 

Case  1.  Adaptive  observer  [6] 

Without  loss  of  generality  let  us  assume  that  the  model  of  the 
dominant  part  of  the  plant  (2.4)  is  in  the  observable  canonical  form 


I  I 


I 

I  o 


x  +  Bu  +  H?i 


•  •  _  i 

yn  *  A^n  +  yA^  B^u 

y  *  c'x  *  [10  ..  0]x  ®  x^  . 


(4.1) 


(4.2) 

(4.3) 


The  algorithm  [6J  for  the  nth  order  adaptive  observer  based  on  the 
dominant  part  (4.1),  (4.3)  without  the  parasitics  (T]  -  0  in  (4. 1) )  is  given  by  the 
equations  (4.4)  through  (4,11),  below.  The  observer  equation  is 


z  ■  Kz  +  [k-a(t)]y  +  b  ( t )  u  +  w  +  r 


(4.4) 


y  ■  c  z  ■  z. 


(4.5) 


where  w  and  r  are  auxiliary  signals  formed  by  the  output  error  e^  ^  y  -  v  and 


the  components 


n-1  ,  ,  n— 2  .  ,  , 

s  +  d0s  +■•••  +  d 
L  n 


Xl’  qi  *  n-l  ,  I  n-2 


(4.6) 


of  the  vectors  v  and  q  as  follows : 

To  1 


v’TA^r 


s  +  d_s  +  . . .  +  d 
2  n 


0 

q'MA2q 


w  ■  -e 


l  v*rAjV 


vTAv 

n 


r  -  -e. 


Matrices  are 


n-j+1 


"dj  "Vi 

•  •  • 

V 

-d 

Cl 

0  -d. 

3 

•  «  • 

♦  f  • 

’dn-l 

0  0 

o  -V 

• 

1  d2 

•  *  • 

Vi 

0  1 

• 

♦ 

CM 

tJ 

dj-2 

• 

•  •  . 

• 

0  0 

.  .  1 

d2 

q'MA^q 

q'MA^q 


_ tl 

0 


-dn  0 


0  o  . 


-Vi 0 


(4.7) 


Tl) 


(4.8) 


}  n-j+1 


and  F  *  7'  >  0,  M  *  M'  >0  while 


K  a  Ik. 


(4.9) 


are  chosen  such  that  c’(sI-K)  d  is  positive  real.  The  adaptive  laws  for 


updating  the  estimated  parameters  are  given  by 
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$  -  -Te^u ■  -a(t)  (4.10) 

*  -  -Me^-  6(t)  (4.11) 

where  <fi  £  a  -  a(t)  and  V  £  b(t)-B  are  the  parameter  errors. 

Case  2.  Adaptive  observer  [5],  [7]. 

The  following  "modal"  canonical  form  is  chosen  for  the  dominant  part 
of  the  plant  (2.4) 


r 1  h'i 

i h 

x  -  a  | -  x  +  Bu  +  H  T| 

I  A 


•  _  ^ 

liT)  -  Afi]  +  uAf  Bfu 

y  »  c'x  ■  x^ 


(4.12) 


(4.13) 

(4.14) 


where  h'»  [11  ...  1],  A  is  an  (n-l)x(n-l)  diagonal  matrix  with  arbitrary  but 
known  constant  and  negative  diagonal  elements  -A^  (i*2...n)  and  a,  B  are  the 
unknown  vectors  to  be  identified.  It  is  shown  in  [5]  that  any  completely  observ¬ 
able  system  can  be  represented  in  this  "modal"  canonical  form.  The  structure 
of  the  adaptive  observer  based  on  (4.12)  with  n  ■  0  is  summarized  in  the  equations 
(4.15)  through  (4.20).  The  adaptive  observer  equation  is 


z  ■  Kz  +  (k-a(t))y  +  6(t)u  +  w  +  r  (4.15) 

y  *  c’z  -  z^  (4.16) 


where  w  and  r  are  auxiliary  signals  formed  by  the  derivatives  of  the  parameter 
error  components  and  the  components 

vi*iTT7xr  qi  "  s"+T"  u  (i=2’-*n)  (4-17) 

i  i 


of  the  vectors  v  and  q  as  follows 
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(4.18) 


Moreover  K 


\  I  h'~ 

:  I  a 

o  i 


-i,  _  i 


is  stable,  the  transfer  function  c'(sI-K)  d  * 


is  strictly  positive  real,  and  d  ■  (10  ...  0)'.  Note  that  the  first  components 
of  the  vectors  v  and  q  are  v^«x^,  =  u*  respectively.  The  adaptive  laws  for 

adjusting  the  parameters  are 


<j>  »  -Fe^u  ”  S(t) 


(4.19) 


f  -  -Me^q  *  b(t) 


(4.20) 


It  is  shown  in  Appendix  I  that  the  stability  of  the  adaptive  observers  in  case 
1  and  2  in  the  presence  of  parasitics  is  equivalent  to  the  stability  of  the 
following  set  of  differential  equations. 


e  ■  Ke  +  d[$'v+  f’q]  -  Hn 


(4.21) 


*i"  *1 


(4.22) 


0  -  -r«jv 


(4.23) 


4*  *  -Mt^q 


(4.24) 


where  K,  d,  v  and  q  are  defined  differently  in  case  1  and  2 

It  is  shown  in  [7]  that  if  u(t)  is  sufficiently  rich  for  an  nth  order 
plant  (ie  it  contains  at  least  n-distinct  frequencies)  without  parasitics 
(Hn *  0)  then  the  system  (4.21)-(4.24)  is  u.a.s.  To  study  the  stability  of  the 


algorithm  with  parasitics  (Hn t  0)  we  express  (4 . 21) — (4 . 24)  in  the  form  of  (2.11) 
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by  introducing  the  composite  error  vector  Z(t)  •  f  where  . 


(4.26) 


The  stability  of  the  homogeneous  part  of  (2.11)  can  be  shown  using  the  same 
Lyapunov  function  as  in  the  case  without  parasitica  [7],  [12]  and  it  is  not 
influenced  by  the  fact  that  A^(t)  depends  on  T|.  However  the  proof  of  u.a.s.  is 
different. 

Theorem  2:  If  u(t)  is  sufficiently  rich  for  the  (n-Hn)th  order  plant  then  the 
homogeneous  part  of  (2.11)  is  u.a.s.  and  the  error  vector  Z(t)  is  bounded  by 
(2.12). 

Proof :  For  u.a.s.  of  the  homogeneous  part  of  (2.11)  a  condition  has  to  be 
imposed  on  u(t) .  This  is  a  consequence  of  the  fact  that  the  components  of  the 
vector  [v'.q']'  have  to  be  linearly  independent  functions  of  time.  This  implies 
that  the  components  of  [x^u]'  have  to  possess  this  independence  property  .  It 

can  be  shown  that  if  u  is  sufficiently  rich  for  the  (n-Hn)-th  plant,  then  the 

components  of  [x,,u>n,],  are  linearly  independent  functions  of  time.  This 
implies  that  the  components  of  [xf,u]'  are  linearly  independent  functions  of 
time  hence  the  homogeneous  part  of  (2.11)  is  u.a.s.  Thus  using  lemma  1  (2.12) 

,  i 

follows . 

Remark  2 ;  As  in  Remark  1  in  this  it  can  also  be  shown  that  u.a.s.  for  the 
homogeneous  part  of  (2.11)  can  be  achieved  for  almost  all  u  which  are  sufficiently 
rich  for  the  nth  order  plant.  That  is  if  u  contains  at  least  n  distinct 

frequencies  except  for  a  particular  combination  for  which  the  condition  of  linear 

independence  of  the  components  of  [xf,u]f  can  be  lost  then  the  homogeneous  part 


of  (2.11)  is  u.a.s. 


NON-MINIHAL  ADAPTIVE  OBSERVER  [4 


[8 


A  non-minimal  state  representation  of  the  plant  (2.4)-(2.6)  is 
(see  Appendix  II  and  Fig.  2) 


where  h'  ■  [1  1  ...  1]  and  A  is  as  defined  in  Section  4. 


Wa^ZICW 


Fig.  2.  Non-minima 1  representation  of  the  dominant  part  of  the  plant. 


L 
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The  structure  of  the  adaptive  observer  for  (5.1),  (5.2)  In  the  absence  of 
parasitics  (le  n  ■  0,  R  » 0)  Is  given  In  [4],  [7]  and  [8]  and  the  basic  equations 

3 

are  reviewed  below.  The  observer  equations  are: 


y  -  a^(fi)y  +  a^Zg  +  b^COu  +  b '  (t)Wg~-  A^y-y) 


(5.4) 


Z  ■  AZ  +  hy 
s  s 


(5.5) 


W  -  AW  +  hu 
s  s 


(5.6) 


where  y(0)  »  0,  Z  (0)  •  0,  W  (0)  »  0.  The  adaptive  laws  for  adjusting  the  un- 
s  s 

known  parameters  are  given  by: 


»  -Te^v 


(5.7) 


¥  -  -Me^p 


(5.8) 


where  <t>  £  [(a  (t)-a-),  (a(t)-a  )’]»  f  »  [  (b_  (t)— b- ) »  (b(t)-b  )']  and 
11  S  1  1  s 

*  /V 

v  *  [y,Zg]',  p  «  [u,Wg]'.  The  stability  of  the  non-minimal  adaptive  observer 
described  by  (5.4)-(5.8)  in  the  presence  of  parasitics  is  equivalent  to  the 
stability  of 


e  -  -X-je  +  d t v ’ 4>  +  q ’ 'i' ]  -  h'R  -h'rTE(O) 


(5.9) 


$  -  -Fe^v 


(5.10) 


¥  •  -Me^p 


(5.11) 


where  (5.9)  is  obtained  by  substracting  (5.4)  from  (5.2). 


Theorem  3:  If  u(t)  is  sufficiently  rich  for  the  (n+m)th  order  plant  then  the 
vector  Z(t)  *  [e^,<fc *  is  bounded,  the  bound  is  of  order  of  u  and  is  given 


bv 


<  uy  —  —  •  A"1! 


lim  I Z(t)l  <  uy  —  —  lA/B-UhllT  Hi  -r^* 

t~  m2  a2  *  £  3  h 


(5.12) 


Proof:*  We  express  (5.9)  through  (5.11)  in  a  compact  form  similar  to  (2.11) 


Z(t)  -  An(t)Z(t)+ARRa(t) 


(5.13) 


where 


An(t) 


I  dv  dp 


-MP  I 


(5.14) 


and  note  that  A^(t)  depends  on  n  through  v. 

The  stability  of  the  homogeneous  part  of  (5.13)  is  proved  in  [4],  [7], 
[8],  [12]  by  choosing  an  appropriate  Lyapunov  function.  However  for  u.a.s.  the 
input  u(t)  has  to  belong  to  a  special  class  of  inputs.  When  n*0  it  is  shown  in 
[7]  that  a  sufficient  condition  for  u.a.s.  of  the  homogeneous  part  of  (5.13)  is 
that  u  has  to  be  sufficiently  "rich"  for  an  nth  order  plant.  Using  the  same 
argument  as  in  the  proof  of  Theorem  2  it  can  be  shown  that  if  u  is  sufficiently 
rich  for  the  (n+m)th  plant  the  u.a.s.  of  the  homogeneous  part  of  (5.1?)  is 
assured.  As  in  Lemma  1  u.a.s.  of  (5.13)  implies  that 


-m.t  t  -m.(t-t) 

IZ(t)l  <  m.e  L  IIZ(0)l  +  f  m.e 
1  0  1 


lA^llR  (r)ldt 


(5.15) 


From  (5.1)  for  some  positive  constants  5^,  we  have 

t  -£2^”^) 

1 R  (t)  II  <  /  C,e  II T  Hit  I  n  (t)  II  dr 

s  o  1  s 


(5.16) 


from  (2.4),  (5.15)  and  (5.16) 


-nut 


ai 


iz(t>i  <v  2  tizco)i  - 

M  4  4 

+  <4 fr ?  VtlARlllR.(0)'  -  ITT01 


•2  .  Z1 

^2"  U 


(5.17) 


-v7“  -V.  lynai  . 

+  •i - ST-4 - 1  - r-  t'n<0»  -  my,  1 

<”2-T>  «2-T> 

®1  £*l  Y 

As  t-*»  lim  II Z  (t)H<uY - —  Ia"  B  JIT  HllA_i.  Note  that  *  A_8  *  Uhl  hence 

t-«o  m2  a2  C2  f  t  s  -R  "R 

(5.12). 


It  Is  pointed  out  that  Remark  2  applies  to  this  case  as  well. 


6.  PARAMETRIZED  OBSERVER  [9] 

It  is  shown  in  Appendix  III  that  the  plant  equations,  (4.1)- (4. 3) 


can  be  written  as 

M(t)  -  FM(t)  +  [Iy,  lu]  M(0)  -  0  (6.1) 

X(t)  -  M(t)p*  +  EXP(F  t)X(O)  +  D(t)  (6.2) 

D(t)  -  FD(t)  +  Hn(t)  D(0)  -  0  (6.3) 

UP  *  Afn(t)  yA^Bfu  (6.4) 

y  *  c'x  -  (6.5) 


where  p *  contains  all  the  unknown  parameters  of  the  matrix  A  and  vector  B  in 
(4.1).  The  adaptive  observer  [9]  for  observing  the  state  x  and  estimating  p* 
when  D(t)  *  0  is  reviewed  below.  The  observer  equation  is 
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x(t)  -  M(t)p(t)  +  EXP(F  t)ic(O) 


(6.6) 


y(t)  -  c'x 


(6.7) 


The  form  of  the  adaptive  lavs  for  updating  the  unknown  vector  p  depend 
on  the  particular  criterion  chosen  for  minimizing  the  output  error  e^  «  y  -  y.  . 
For  the  first  adaptation  scheme  the  error  criterion  is 


?  ,Ie2 
41  2  1 


(6.8) 


and  the  adaptive  law  [9],  [10]  is  given  by 


p(t)  *  -GM' (t)c  e. 


(6.9) 


where 


G  *  G  >  0. 


For  the  second  adaptive  scheme  the  error  criterion  is 


"  /  [c'M(T)p(t)  +  c'EXP[Ft]x(0)  -  y(t) ]e"^ ^t-T^dT 


(6.10) 


where  q  is  a  positive  constant,  and  the  adaptive  laws  [9]  are 


p  *  -G [R(t)p(t)  +  r(t) ] 


(6.11) 


R(t)  -  -qR(t)  +  M' (t)cc 'M(t) 


R(0)  -  0 


(6.12) 


t(t)  -  -qr(t)  +  M'(t)c[c'EXP[Ft]x(0)  -  y (t) ]  r(0)  »  0  (6.13) 


We  now  analyze  the  stability  of  each  of  the  two  adaptation  schemes  in  the 
presence  of  the  parasitic  input.  We  also  derive  bounds  for  the  composite  identi¬ 
fication  and  observation  error 


x(t)  -  x ( t ) 


p(t)  -  p* 


(6.14) 


In  the  first  adaptation  scheme  the  state  error  e  and  parameter 


error  Ap  satisfy 


e  -  -M(t)GM'  (t)cc'e  +  [FH(t)  +  (Iy,Iu)]Ap  +  exp(Ft)Fe(0)  -  FD(t)-Hn(t) 

Ap  “  -GM'  (t)cc'e 


(6.1 

(6.1 


Combining  (6.17)  and  (6.18)  the  equation  for  the  composite  error  Z(c)  is 


Z (t)  -  A  (t)Z(t)  +  B  (t) 
n  h 


(6.1 


where 


An(t) 


-M(t)GM' (t)cc’  FM(t)  +  [Iy,Iu] 


-GM' (t)cc’ 


(6.1 


and 


Bn(t) 


Ft, 


e  Fe(0)  -  HTl(t)  -  FD(t) 
0 


(6.1 


The  following  Theorem  establishes  the  condition  for  the  u.a.s.  of 
the  homogeneous  part  of  (6.17)  and  gives  a  bound  for  Z(t). 

Theorem  4:  If  u(t)  is  sufficiently  rich  for  an  (n+m)th  order  plant  then  the 
nth  order  adaptive  observer  given  by  (6.6),  (6.7),  (6.9)  is, stable  in  the 
presence  of  the  parasitic  part  (6.4)  of  the  plant  (6.1)— (6.5)  in  the  sense  that 
the  composite  error  vector  Z  is  bounded.  A  bound  on  Z  as 
t-*co  is  of  order  of  y  and  is  given  by 

f. 


lim  «Zl  <  uy  —  —  H  Hit  II  A~1B  [1  +  8  Fll] 

t-  «2  m2  C  f  f2 


(6.2l 


Proof :  If  u  is  sufficiently  rich  for  the  (n+m)th  plant  then  the  components  of 
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M'(t)c  are  linearly  independent  functions  of  time  [9].  This  implies  that 
there  exist  constants  k^,  k2  and  T  such  that 


t+T 

0  <1^1  <  ft  M'  (t)cc  'M(x)dx  <  K_I  for  all  t  >  0 


(6.21) 


is  satisfied.  Using  the  same  proof  as  in  Theorem  1  of  [9.]  it  can  be  shown  that 
if  (6.21)  is  satisfied  then  the  homogeneous  part  of  (6.17)  is  u.a.s.  with  the 
rate  of  convergence  no  less  than  m2  where 


K  min  X [G] 

l2  "  - - T 

[1  +  nK  max  X[G]r 


Thus  from  (6.17) 


(6.22) 


”m2t  »■  -nu(t-T) 

D Z (t) I  <  m_  «.  HZ (0)1  +  fm.  £  Hb  (T)fldr 

o  l  n 


(6.23) 


Bn<c) i  ^ flA  2  llFe(0)l1  +  a Hil U n (t) II  +  lFllD(t)l 


t  -f,(t-t) 

»D(t)ll  <  /‘f  a  2 
o  1 


II  Hll  II  n  (t)  II  dr 


(6.24) 


(6.25) 


From  (6.23),  (6.24),  (6.25)  and  (3.4) 


-nut  F2C  _m2C 

I Z ( t ) H  <  m.  i  2  I Z (0) I  +  p,m  — — - 2 

1  1  1  (m2  -  f  ) 


+  P2m2 


-a7t/u  -m_t 

(A  l  ) 


<"2  - 


-m0t 

(1-*  2  ) 
+  P3ml  —  2' 


(6.26) 


where  £±  are  positive  constants,  f2  <min  |x[F]|  and  P]_,  p2>  are  given 


P,  *  (f,«Fe(0)l  +uy  —  1  AgXB  J  f  ( - 

1  1  cu  f  f  1 v  a 


*  -  “r4  -  ^  1  «-2n 

VT  2  <f2  - 


P2  [Vl a -  +  \  'HI  <1(0)1  -  UY  -I  lA  *B  I  (  - i_  +  lHl)]  (6.28) 

(f  2  3 2 

U2  u  '  f 2-“ 


p,  -  UY  —  JA,  B,UhI  [1  +  »F»] 

3  o2  f  f  f2 

As  t-*®  (6.20)  follows. 


In  the  second  adaptation  scheme  the  state  and  parameter  errors 


satisfy 


e  -  [FM(t)  +  (Iy,Iu)-  M(t)GR(t)  ]^p  -  M(t)G  zV  (t)cc  ' V  re(0)£ 

o 

+  M(t)G  zV  (T)cc,D(T)rq<t_T)dr  +  JlFtFe(0)  -  FD(t)  -  Hn(t) 
o 


dt  +  r  Fe(0)  -  FD(t)  -  Hn(t) 


Ap  -  -GR(t)Ap  -  G  ZCM’ (T)cc,lFTe(0)l"q(,:'T)dT  +  G  ZCM'  (r)cc  'D(T)jTq(t’T)dT 

°  0  (6.3 

respectively. 

By  defining  the  composite  error  Z  as  in  (6.14)  it  can  be  easily 
shown  that  the  stability  of  the  second  adaptation  scheme  is  equivalent  to  the 
stability  of 


Z  »  A  (t)Z  +  E(t)  +  C(t)  +  K(t) 

n 


(6.32: 


where 


An<c>  - 


(FM(t)  +  (Iy,Iu)  -M(t)GR(t)] 


-  GR(t) 


lFtFe(0)  -  M(t)G  ZC M'  (T)cc'«.FTe(0)  l‘q(t"T) 
o 

-  G ZCM' (T)cc'ZFTe(0)i~q(t_T) 


-  FD(t)  -  H  T;(  t )  +M(t)GZt  M’  (r)  cc  ’  D(t)  rq(t~T )  dt 

o 


(6.33: 


(6.34: 


(6.35: 


G  ZC M*  (t)cc’D(t)  rq<‘t_t)dr 


(6.36 
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Note  that  E(t)  vanishes  with  an  exponential  rate  not  less  than  min  {f 2 * q  J  ■  w^. 

The  u.a.s.  of  the  homogeneous  part  of  (6.32)  and  a  bound  on  the 

composite  error  Z(t)  are  established  by  the  following  theorem. 

Theorem  5:  If  the  input  u  is  sufficiently  rich  for  the  (n+m)th  order  plant 

then  the  nth  order  adaptive  observer  specified  by  (6.6),  (6.7)  and  (6.11) 

through  (6.13)  is  stable  in  the  presence  of  the  parasitic  part  of  the  plant. 

The  bound  on  Z(t)  as  t-**>  is  of  order  of  y  and  is  given  by 

f  f3 

lim  HZ(t)8  <yy  — BA:1BJ8HI[1+Ti  IfU  +4  —  [U+T1  *Bl)6 
t~  m2  a2  f  f  f2  f*  q  C2 

2 

®i  ?■]  2  II  rll  ^1  £1  I  _i 

+  PY  »A  1BJ»H«]2+  -^-^-[(l+Ti  BBI)6  +  PY  7^' 8  A /b  J  1  H«  ]  ]  (6.37) 

a2  h  f  f  f^  q  h  a2  h  f  f 

Proof:  Since  u(t)  is  sufficiently  rich  for  the  (n+m)th  order  plant  the 
components  of  M'(t)c  are  linearly  independent  [9].  This  condition  guarantees 
the  existence  of  constants  K  and  T  such  that 


(T)cc'M(t)d->KI  >  0 


for  all  t>0 


(6.38) 


In  Theorem  3  of  [9]  it  is  shown  that  if  (6.38)  is  satisfied  then  the 
homogeneous  part  of  (6.32)  is  exponentially  stable  with  a  rate  of  convergence 
*  mintq.t^l  where 


02  *  K2  min  A [ G ] 


(6.39) 


Thus  for  some  positive  constants  w^,  s^,  G^,  5  and  it  can  be 


shown  that 


I  Z  ( t )  II  <  m-e  2  IIZ(0)H  +  m,  /  l  [ »  E  ( t )  H  +  II  C  (  t  )  II  +  I  K(t)  H  ]dr(6,40) 

1  1  0 


I  E ( t )  II  <  wxi 


(6.41) 


IC(t)l<s.l  ‘  +  MY  —  lA/’B  J I  Hi  [1+7^-  1  Fl  ]  +  lM(t)llK(t)l 

x  <*2  t  t  1 2 

^1  ®i  £-i  _i 

lM(e)l  <  c.A  1  +  [(1+-==-  lBl)6  +UY  —  ■=*■  iHl I A-^B-I  ] 

1  f2  h  a2  h  f  f 

_62C  fl  °1  I  cl  -l  €1 

«K(t)i  <  6.e  2  +  py  -y  —  —  Ia/b,HhI  [(1  +  7^-  lBl)6 

1  f2  °2  q  f  f  *2 

+  UY  iil  IHIIA^BJ] 

“2  52  f  f 


(6.42) 


(6.43) 


(6.44) 


Then  from  (6.40)  through  (6.44) 

2 

—n*}t  m.  f  i  H  r~> | 

I  2(t)l  <  n  2,  2  +  UY  —  •a"1B,I«hI  [1  + IF* +-i  —  ((1+7^-  iBl)6 

•  1  m2  2  f  f  f2  f2  q  ^2 

3  2 

al  D  rj  ai  _i  2 

+  UY  —  t*-  II A  XB  I II  HU  )  +  —  ((1+7=-  lB*)6  +  uy  ~  7“  *a/b  JlHl)2]  (6.45) 

^2  ^2  ^  Q  ^2  a2  ^2  ^  ^ 

.  02  «2 

where  s2*min[—  , f 2 ] ,  e2»min[—  ,C2,f2] 

62  -  min[e2,q],  n2  -  min[s2,e2,62] , 


C2<  min|X[A]|  and  6  *  sup  u(t)  for  all  t  >  0.  As  t-*-”  (6.45)  reduces  to  (6.37) 

Remark  3:  In  both  first  and  second  adaptation  schemes  u  was  required  to  be 
sufficiently  rich  for  an  (n-Hn)th  order  plant  in  order  for  conditions  (6.21)  and 
(6.38)  to  be  satisfied.  However,  as  it  was  pointed  out  in  Remark  2,  for  almost 
all  u  sufficiently  rich  for  an  nth  order  plant  the  results,  obtained  will  still 
be  valid. 


7.  DISCUSSION 


In  all  the  error  bounds  derived  In  the  preceeding  sections  the 
common  factor  is 


fv  ~  rp  IIaI^JIIIhII  (7.1) 

“2  *2  f  f 

The  most  important  terms  appearing  in  (7.1)  are  k*,  the  "speed  ratio"  of  slow 
vs  fast  phenomena,  and  Y  which  is  a  characteristic  of  the  input.  The  dependence 
ot  the  error  bound  on  H  shows  that  the  adaptive  schemes  considered  are  robust 
with  respect  to  the  parasitics  in  the  sense  that  as  H'-’O  the  error  bound  goes 
to  zero.  We  demonstrate  the  effect  of  H>  and  Y  by  digital  simulation  of  the 
adaptive  observer  (Section  4,  Case  1)  for  the  plant. 


r  h 

—  — < 

— * 

-5  1 

X  + 

0.9 

Xf  + 

1.45 

-10  0 

0.5 

1 

2.25 

—  _ 

—  — 

Rx^  ■  -4x^  “  2u 
y  -  [1  0]x 


(7.2) 


(7.3) 

(7.4) 


become 


Using  the  transformation  ^lax^  +  0.3u  the  plant  state  equations 


(J.T]  »  -4?,  +  O.SF-u 


(7.5) 

(7.6) 


y  -  [1  0]x 


(7.7) 


The  adaptive  observer  tor  (7.5),  (7.7)  with  11  ■  0  is 
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-6 


-8 


i|  U.vtn  fb,^)  -o  o 

*  +  ^  y+Ls  u”ei  •  *i 

°J  L®"a2(t)J  Vfc)  vrV 


where 


A. 


the  components  of  the  signals  v  and  q  are  generated  by 


v2  -  -3v2  +  -  Vj^,  -  -3q2  +  u  -  qt 


The  adaptive  laws  for  adjusting  the  parameters  are  taken  as 


a(t)  -  TejV,  b(t)  -  -Me^q 


(7.8) 


P  -1.  r-p“  "1».F  "1 

L°  lJ  L°  75J  L°  7-8J 


and 


(7.9) 


(7.10) 


The  dependence  of  the  error  bound  on  H*  is  illustrated  in  Figs.  3  and  4.  For 
^*0.2  and  input  u»5sint  +  5sin  2.5t  that  is  Y  »18.2,  the  observation  error  e ^ 
is  relatively  small.  However  the  parameter  errors  are  significant:  10.4%  tor 
a^(t)  and  12%  for  b^(t).  Reduction  of  ^  by  a  factor  4  that  is  **»0.05  results 
in  a  reduction  of  the  parameter  errors  by  approximately  the  same  factor  as 
shown  in  Fig.  4b,  c.  The  observation  error  e2  is  almost  zero  in  this  case 
(Fig.  4a). 

To  examine  the  effect  of  Y  «suplu(t)|  on  the  error  bound  the  value 
of  is  kept  the  same  as  in  Fig.  4  but  the  input  is  changed  to  u  »  5sint  + 

+  15sin  2.5t  that  is  Y  is  increased  to  Y  -42.5.  The  results  obtained  are 
shown  in  Fig.  5.  It  is  clear  that  the  value  of  Y  is  crucial  for  the  error  bound. 
Increasing  Y  by  a  factor  of  2.3  results  in  an  increase  of  the  parameter  error 
by  a  factor  of  about  10.  Moreover  the  observation  error  although  bounded,  is 
oscillatory  and  not  close  to  zero. 
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(a)  (b) 

Fig.  3.  u  -  0.2,  a  -  5  sint +5  sin  2.5t  (y-18.2) 
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Fig.  4.  u  *  0.05,  u  *  5  sint+5  sin  2.5t  (y*18.2) 
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Fie.  5.  u*G.05.u*5,  sint+  15  sin  2.5t  (y  *  42.5) 


The  above  simulation  results  show  that  the  choice  ot  the  excitation 
Input  signal  is  critical  in  adaptive  schemes  with  modeling  errors.  The  effect 
of  the  input  on  the  error  can  be  even  more  crucial  than  the  effect  of  K 
Rapidly  varying  inputs  (high  Y)  will  result  into  bad  estimates  of  parameters 
and  states.  Inputs  of  this  class  excite  the  parasitic  part  of  the  plant 
considerably  and  have  adverse  effects  on  the  results  of  adaptive  schemes  with 
modeling  error.  Caution  is  needed  in  selecting  an  input  excitation  signal. 

Our  results  show  that  the  richness  condition,  should  be  satisfied  but  with  an  in 
put  which  has  a  low  value  of  Y.  Slowly  varying  inputs  are  appropriate  as  long 
as  they  do  not  reduce  the  convergence  rate  considerably. 

Convergence  rate  with  which  Z(t)  exponentially  decays  in  Z -A^,(t)Z 

is  sensitive  with  respect  to  Y  only  for  very  low  values  <*f  Y .  Thus  we  can 

improve  the  error  bound  by  keeping  Y  as  low  as  possible  for  not  to  be 

affected  appreciably.  From  (7.1)  it  is  clear  that  larger  m^  reduces  the  error 

bound.  However  depends  on  the  input  u  and  on  the  adaptive  gains  and,  hence, 

for  low  values  of  Y  there  is  a  trade-off  between  Y  and  m^,  that  is  m^  cannot 

be  improved  through  the  choice  of  the  input.  Its  improvement  by  the  choice 

of  the  adaptive  gains  can  only  be  done  by  trial  and  error  since  in  all  the 

adaptive  schemes  considered,  except  in  the  parametrized  observer,  m^  is  not 

explicitely  related  to  the  adaptive  gains.  In  the  case  of  the  parametrized 

adaptive  observer  the  expression  for  gives  more  information  about  the 

dependence  of  the  error  bound  on  other  quantities.  For  the  first  adaptation 
kjminX[G] 

scheme  nu  ■  - r  .  In  this  case  the  best  we  can  do  to  improve 

[1  +  n  k2maxX[G]]^ 

m^  and  consequently  the  error  bound  is  to  make  maxX[G]  *  min  X[G] .  The 
dependence  of  on  the  order  n  of  the  dominant  part  of  the  plant  indicates 
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that  the  error  bound  will  be  higher  for  a  higher  order  dominant  part  of  the 

-qT 

plant.  For  the  second  adaptive  scheme  ■minfq^]  where  02  «Ke  min  a.[G]  . 

In  this  case  can  be  increased  arbitrarily  by  increasing  G  and  choosing  q 

appropriately.  However  the  bound  for  this  scheme  is  proportional  to  IIgII  and 

it  is  not  clear  whether  an  increase  of  m2  through  G  will  reduce  the  error. 

The  term  m,  in  (7.1)  indicates  the  dependence  of  the  bound  on  the  Initial  error 

a 

vector  Z(0).  The  factor  IIa^b^II  in  (7.1)  depends  on  the  characteristics  of 
the  parasitic  part  and  the  term  HhII  is  a  measure  of  the  coupling  between  the 
dominant  and  parasitic  part  of  the  plant.  Apart  from  the  common  factor  (7.1) 
the  bounds  obtained  for  the  parametrized  adaptive  observer  are  also  functions 
of  the  characteristics  of  the  observer  gain  F.  The  bound  for  the  second 
scheme  is  more  complicated  due  to  the  complexity  of  that  scheme. 

8.  CONCLUSIONS 

The  results  of  the  paper  show  that  if  the  plant  is  stable,  the  adapta- 
tive  schemes  considered  remain  bounded  despite  the  reduction  of  the  model  order. 

The  bound  on  the  observation  and  parameter  error  is  of  the  order  of  the  singular 
perturbation  parameter  ^,and  is  also  a  function  of  the  characteristics  of  the 
input,  the  initial  parameter  and  state  error  and  the  convergence  rate  which 
would  be  achieved  if  there  were  no  modeling  error.  The  dependence  of  the  error 
on  the  input  characteristics  is  found  to  be  crucial  and  the  most  desirable 
excitation  signals  are  those  which  are  sufficiently  rich,  but  have  a  low  value 
of  Y  *sup|u(t)i.  A  trade  off  between  Y  and  m2  should  be  made  when  selecting  the  input 
sigt  il.  The  input  does  not  have  to  be  sufficiently  rich  for  the  full  order 
plant.  The  results  of  the  paper  are  valid  for  almost  every  u(t)  suffici«",'1y 
rich  for  the  nth  order  dominant  part  of  the  plant.  Thus  u(t)  can  contain  as 
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many  frequencies  as  required  Co  be  rich  for  Che  n -order  plane  excepC  for  a 
particular  combinacion  of  frequencies  for  which  Che  richness  of  u(C)  is  reduced 
by  Che  parasiCic  input  7](t).  Extensions  of  these  results  Co  model  reference 
adaptive  control  and  ocher  closed  loop  adaptive  schemes  is  a  topic  for  future 
research. 
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APPENDIX  I 


The  stability  of  the  adaptive  observer  in  the  presence  of  parasitics 
for  Case  1  and  2  is  equivalent  to  the  stability  of  the  following  set  of 
differential  equations 


e  ■  Ke  +  $x^  +  Fu  -  Hn  +  w  +  r,  e^  ■  c'e 


(1.1) 


$  *  Te^u 


F  ”  -Me^q 


(1.2) 


(1.3) 


Proposition;  For  some  vector  signals  v,  q,  w  and  r  with  v  ■  G(p)x^,  q  *  G(p)u, 
•  • 

w  •  w($,u)  and  r  ■  r(F,q)  the  system  (1.1)  is  input  ( [x^,u, n ' ] ' )-output  (e^) 
equivalent  with  thesysteml.4  provided  (c',K)  is  completely  observable 


e  »  Ke  +  d[4>'v  +  F'q]  -  Hn,  e^  -  c'e  *  e^ 


(1.4) 


Proof ;  From  (1.1),  (1.4) 


c'(pI-K)  X[$x^  +  Fu  +  w  +  r  -  d(4>'v  +  F'q)] 


(1.5) 


where  'p'  is  the  d/dt  operator.  From  (1.5)  we  have 


Z  pn-± [ <>ix1  +  ^  +  wi  +  ri  “  di<*),v  “  djF'qJ  »  0 


(1.6) 


where  i  denotes  the  ith  element  of  the  corresponding  vector.  (1.6)  is  satisfied 
by  choosing  v,  q,  w  and  r  as  given  in  Case  1  and  Case  2.  By  considering  the  two 
equivalent  systems  (1.1),  (1.2),  (1.3)  and  (1.4),  (1.2),  (1.3)  .  it  is  obvious 
that  boundeness  of  [s' ,F* ] '  will  imply  boundeness  of  [e' ,F' ] . 
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APPENDIX  II 


From  (5.1)  and  (5.3) 


5$  ■  °<*>  ■  y>> + 


where 


G  (s)  -  c'(sI-A)  B 
P 


is  the  transfer  function  of  the  plant  when  Hn  ■  0.  Let 


bn-1  *  i  n-2  ,  ,  . 

« s  +  b0s  +  •  #  •  b 

p 


-  ,  n-l  ,  n-2  . 

s  +  a,s  +  ous  +  ...  +  a 
12  n 


where 


“  “  C°i  »a2’  *  *  *  »anl  *  and  B  *  Cb1.t>2,. 


are  the  unknown  parameter  vectors. 


Consider  a  polynomial  tt  (s+X1)  which  is  relatively  prime  to  the 

i-2 


numerator  as  well  as  the  denominator  polynomials  of  G  (s)  and  A.  j4  AJ  for  i, 

P  1  i 

j*2,3,...,n).  Dividing  the  numerator  and  denominator  of  G(s)  by  it  (s  +  A  ) 

i=2  i 


and  expanding  them  into  partial  fractions  we  have 


-1, 


c’adHsI-A)  H  n(s) 


bl  +  s  +  A, 


+  *  .  + 


s  +  A 


G(s) 


^2(3  + Ai) 


u(s) 


(II 


s  -  a„  - 


1  s  +  A, 


s  +  A 


s  -  a„  - 


Note  chat 


n 


1  s  +  *  *  *  s  +  A 

2  n 


-1 


c  Tadj (si- A)  x  Is 


n-l  n-2 


ZJLULlU 


tZ2(s  +  x±) 


i~2(s  +  Xi) 


iA 


n  tli 


‘2i 


i*2  s  +  A’  i*2  s  +  A. 

i  1 


5  Cni 


i*2  s  +  A . 


(II 
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II. 1  can  be  written  as 


1  n  H>,u(s)  +  *<y(s)  +  T.HTl(s)] 

y(s)  -  7[b1u(s)  +  a^y(s)  +  ±£2— - ;TT - 


(II. 3) 


where 


*CU*  t2i,"’,tni^ 


II. 3  is  represented  In  a  block  diagram  form  in  Figure  2.  The  block  diagram  of 
Figure  2  contains  (3n-2)  integrators  and  is  a  nonmlnimal  realization  of  the 
dominant  part  of  the  plant. 

The  term  h'exp(A  t)x(O)  in  Figure  2  is  added  so  that  Figure  2  is 

equivalent  to  the  corresponding  figure  of  a  minimum  realization  of  II. 3  including 

initial  conditions.  Here  x  »  [x0 ,x_ , . . . ,x  ]  and  x  is  the  state  of  the  minimal 

•  l  j  n 

realization  based. on  II. 3. 

The  nonmlnimal  state-space  representation  given  by  (6.1),  (6.2)  can 
be  easily  obtained  from  Figure  2  by  defining 


Rg  -  [r2,r3, . . . ,rQ] ' ,  wg  -  [w2,w3, . . . ,wnJ ' 
'Zg  -  [z^Zy  . . .  ,zn]  ’  and  Tg  *  [T2,Tj,  . . .  ,T^1 ' 


I. 


I 

L 
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APPENDIX  III 


where 


Equation  (2.4)  can  be  represented  as 


i  ■  Fx  +  gy  +  Bu  +  Hn 


-fx  1  0  ...  0 

-f2  0  1 

:  :  i 

•  • 

-f  0  0  ...  0 


l 

and  g  satisfies  gc'  *  A-F. 
From  III.l 


x  -  *Ftxo  +  aF(t_T)[gy(T)  +  Bu(x)]dT  +  fl  *F(t"T)Hn(T)dr 


The  first  convolution  integral  can  be  reduced  to 


fC  £F(t-T)[Iy(T)>  Iu(T)jdT.  p*  .  M(t)p* 


where  p*'  -  [g',B'].  Thus  (6.1)-(6.3)  follows  by  taking 


D(t)  -  /C  JiF(t'T)Hn(T)dT 


(III.l) 


(III. 2) 


(III. 3) 


(III. 4) 


(III. 5) 


i 
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